Microscopic theory of the inverse Faraday effect 



in 
o 
o 

(N 
D 



o 
o 



> 
o 
\o 
o 

ON 

o 
vn 
o 



o 
o 



X 
J3 



Riccardo HerteQ 

Institute of Solid State Research (IFF), Research Center Jiilich, D-52425 Jiilich, Germany 

An analytic expression is given for the inverse Faraday effect, i.e. for the magnetization occurring 
in a transparent medium exposed to a circularly polarized high-frequency electromagnetic wave. 
Using a microscopic approach the magnetization of the medium due to the inverse Faraday effect 
is identified as the result of microscopic solenoidal currents generated by the electromagnetic wave. 
In contrast to the better known phenomenological derivation, the microscopic treatment provides 
important information on the frequency dependence of the inverse Faraday effect. 



I. INTRODUCTION 

While the equations describing the transfer of energy 
and momentum to a medium by an electromagnetic wave 
are commonly known from undergraduate students' tex- 
books, the study of the transfer of angular momentum 
of an electromagnetic wave has received relatively little 
attention. The first prediction that circularly polarized 
light should cause an angular momentum-flux is due to 
Sadovskii and dates back to 1897. The original work of 
Sadovskii is written in Russian and is difficult to find 
(see footnote 1 in Ref. 0). The first detection of the 
Sadovskii effect, i.e. the mechanical measurement of a 
torque exerted by circularly polarized light on a half- 
wave plate, has been reported in 1936 by Beth l] and 
Holbourn independently. The interesting history of 
the study of angular momentum of electromagnetic waves 
is reported in several review papers Since magnetic 

moments have an angular momentum, the possibility to 
influence the magnetization in a magnetically ordered 
material by means of circularly polarized light appears 
to be a straightforward continuation of the old experi- 
ments by Beth and Holbourn. It is therefore remarkable 
that only very recently, about 70 years after these his- 
toric experiments, it has been demonstrated by Kimel et 
al. [3 that spins can indeed be manipulated with a circu- 
larly polarized laser beam. Apart from their importance 
for fundamental processes concerning the magnetization 
dynamics, these spectacular experiments are also inter- 
esting, because of their potential technical implications. 
These studies may pave the way towards an ultrafast, 
focussed laser-controlled magnetic writing process that 
could eventually replace the comparatively clumsy data 
storage technique being in use nowadays, which involves 
microscopic coils and write heads that may generate un- 
desired fringing fields. The effect reported by Kimel et al. 
is attributed to the generation of a stationary magnetic 
field resulting from irradiation with circularly polarized 
light, which is known as the inverse Faraday effect. 

The inverse Faraday effect (IFF) has first been pre- 
dicted by Pitaevskii more than 40 years ago 0, while 
the name (IFF) seems to have been coined only four years 



later by van der Ziel et al. 0, Q . The IFF was originally 
derived on a phenomenological basis and has received a 
certain publicity among physicists, because it was briefly 
treated in some older editions of a well-known textbook 
series on theoretical physics 0]. A much more detailed 
description of the IFF has been elaborated in the late sev- 
enties and the early eighties by theoretical plasma phvsi- 
cists, mainly in the former Soviet Union pJl IllL Il3. Il3| . 
This literature is presumably hardly known in the part 
of the scientific community that is currently investigat- 
ing ultrafast magnetization processes. The scope of this 
paper is to provide a simple derivation of the IFF that in- 
volves microscopic currents, described in more detail by 
Karpman et al. 'l3l, rather than a purely phenomeno- 
logical approach. In the wake of the recent studies by 
Kimel et al., this short review of the theory of the IFF 
is likely to be of interest for researchers working in the 
domain of magnetism who may not happen to be experts 
in theoretical plasma physics. 



II. BASIC EQUATIONS 

When a material is irradiated with a high-frequency 
(HF) electromagnetic wave, the primary coupling be- 
tween the wave and the material is given by the inter- 
action between the electrons and the electric field of the 
wave. One may therefore neglect in a first approximation 
the wave's fluctuating magnetic field SB. 

For simplicity, we also neglect any quantum mechan- 
ical effects and assume that the band structure of the 
material is unimportant for the motion of the electrons. 
In the case of a metallic material the conducting elec- 
trons can be treated as a collisionless plasma in which 
the electrons can move freely, at least on the time scale 
given by the period of the HF field. Such an electron 
plasma can conveniently be represented as a fluid with 
density n{r,t) and velocity v{r,t). The electron density 
n and the velocity v{r,t) are related to each other by the 
continuity equation 
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V {nv) = 



(1) 



The electron density n and the velocity v give rise to an 
electric current density j according to 



'Electronic address: "r . her t el@f z- j uelich. de| 



(2) 



2 



which, in the simplest approximation, is proportional to 
an electric field E 



j = aE. 

In the case of an oscillating electric field 
SE{x, t) — E exp{ikr — iujt) 



(3) 



(4) 



of high frequency lu and wave vector fe, the conductivity 
a of an isotropic, coUisionless plasma is 



(5) 



where m is the electron mass. The brackets ( ) denote a 
time average of the electron density n over several periods 
of the HF field. 

For the following, only the temporal oscillation of the 
local field at a given point r will be required, so that the 
spatial dependence of the electromagnetic wave SE can 
be omitted. Note that the amplitude E is in general a 
complex quantity. The physically relevant real value of 
the field is given by 

K {SE} = ^E exp{-iujt) + c.c. , (6) 

where c.c. denotes complex conjugation of the previous 
term. 



III. SEPARATION OF TIME SCALES 

To investigate the influence of an HF field on a plasma, 
it is useful to treat different time scales separately. The 
time scales are split by representing the electron density 
n, the velocity t), and any other quantity a which is re- 
quired for the description of the plasma as a sum of two 
parts 



a = (a) + Sa 



(7) 



one of which ((a)) is the time-averaged value that is con- 
stant or slowly changing on the time scale given by the 
period of the HF field, and the other {6a) is an oscillating 
part with the same time dependence oc exp(— iwt) as the 
HF field HI. 

We may safely assume that the magnitude of the os- 
cillating part 5n of the electron density, i.e. the size of 
the density fluctuations induced by the HF field, is much 
smaller than the stationary value of n: 



Sn <C {n) 



(8) 



Inserting the separation of stationary and oscillatory 
quantities into the continuity equation ^ yields 



d_ 

dt 



{{n) + (5n) + V i{n){v) + {n)5v + 5n{v) + SnSv) = 0. 

(9) 



Owing to {5a) = 0, the time-averaged continuity equation 
is 

^^{n) + V{{ji){v) + {6nSv))=0 . (10) 

This is the so-called slow time scale of the continuity 
equation. Using the approximation of eq. 0, the fast 
time scale reads 



d 

—5n + V {{n)6v) = 
at 



(11) 



The same approximation yields the oscillating part 5j of 
the current density j as 



5j — e{n)Sv 
Hence, the fluctuating part of the velocity is 



Sv = 



[n)e 



SE 



(12) 



(13) 



According to eq. (|ll|l and eq. H13|) . the fluctuating part 
of the electron density n is given by 

dn = — ^ V {aSE) . (14) 

The slow time scale term of the current density (j) is 

(j) ^ e{{n){v) + {SnSv)) (15) 

or simply 

(j) = e{Sn5v) (16) 
if (v) — 0, which shall be assumed henceforth. 

IV. INVERSE FARADAY EFFECT 

Equation (|16|l describes the occurrence of a station- 
ary current density {j) that is solely due to the high- 
frequency oscillations of the electron density and veloc- 
ity induced by the HF field. This term will eventually 
lead to an expression for the IFE, i.e. for the magneti- 
zation of the plasma in the field of a circularly polarized 
electromagnetic wave. 

To further evaluate the right-hand side of eq. ifTCI) it 
should be beared in mind that the term Sn ■ dj describes 
the product of two real quantities. Therefore, 

Sn ■ Sv — 5ft {nexp(— icji)} ■ |j cxp(— i[x)i)| 

_ 1 
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nv* + n*v + nj exp(— 2icji) + n*v* exp{2iLut) 

(17) 



and hence 
{SnSv) 



— {nv* + n*v) 
4 



4e^(n}tj 



'E* -V aE]- 



(18) 
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Using the identity 

A-VB-B-VA = Vx{AxB)- (BV) A + (AV) B 

(19) 

one obtains 



a*E* X aE 



with 



4:e{n)uj 



r (20) 



(21) 



The term T describes the so-cahed ponderomotive force 
acting in the plasma. The ponderomotive force leads to 
a special type of stationary currents generated by a HF 
field. These currents result from spatial inhomogeneitics 
of the plasma or of the field of the wave. This term F 
is discussed in detail, e.g., in Ref. Although it can 
be assumed that ponderomotive forces also play a role in 
the recent experiments by Kimel et ai, this term does 
not describe the IFE. It is the first term on the right- 
hand side of eq. H2U|I that has the form of a solenoidal 
magnetization current j^: 



(22) 



Hence, the magnetization M generated in the plasma by 
the HF field is 



M = 



2 



a*E* X aE 



I ecu" 



ExE* 



(23) 
(24) 



where Wp = (47r(n)e^/m)^/^ is the plasma frequency. 



V. DISCUSSION 

It should be noted that eq. (|24() has been derived with- 
out making any assumption concerning the polarization 
of the electromagnetic wave. In the case of a circularly 
polarized wave propagating in the z direction, the last 
term on the right hand side can be written as 



ExE* ^ ±i\E\ 



(25) 



while it is equal to zero in the case of linear polarization. 
The plus and the minus sign in eq. (|25|l refers to left and 
right circular polarization, respectively. Therefore, M is 
a real-valued, stationary magnetization that is induced 
in the medium by a circularly polarized electromagnetic 
wave. The magnetization is parallel to the axis of prop- 
agation of the wave and its sign depends on the chirality 
of the wave. 

If the IFE is applied to change the magnetization in 
magnetically ordered materials, the magnetization in- 
duced by the HF field according to eq. H24|l should re- 
orient the magnetic moments in the sample, leading to 



a permanent magnetization that persists after the HF 
field is switched off. The mechanism that converts the 
field-induced magnetization currents into a permanent 
magnetization is not yet clear. The simplest explanation 
would be the magnetostatic alignment of the magnetic 
moments in the field generated by the IFE. However, a 
more complicated mechanism on an electronic level can- 
not be ruled out currently. 

In the derivation of the IFE according to eq. (|24|l it 
has been assumed that the electromagnetic wave is not 
absorbed by the medium. Since the value of the high- 
frequency conductivity is purely imaginary, cf. eq. lO, 
the induced current density 5j is phase-shifted by 7r/2 
with respect to the electromagnetic wave. Hence, the 
power density transferred to the plasma 



{6j6E)='^\E\' 



(26) 



is also purely imaginary, which means that the electro- 
magnetic wave is not absorbed. However, it should not 
be concluded from this absorption-free derivation of the 
IFE that the absorption of a circularly polarized electro- 
magnetic wave does not lead to a magnetization of the 
medium. 

Interestingly, none of the historical arguments on the 
transfer of angular momentum mentioned in the intro- 
duction has been used for the derivation of the IFE. In 
fact, the investigation of the angular momentum trans- 
ported by a circularly polarized electromagnetic plane 
wave leads to several complications. It has been demon- 
strated, that an infinitely extended, circularly polarized 
plane wave, surprisingly, does not carry any angular mo- 
mentum . This counterintuitive result is not in contra- 
diction to the experimental proof of the angular momen- 
tum demonstrated by Beth Q and Holbourn , since it 
has also been shown that a real, coUimated beam of cir- 
cularly polarized light does have an angular momentum 
flux [31 . Using higher order terms in the so-called geo- 
metric optics approximation 15], an analytic expression 
for the spin fiux density of a coUimated circularly polar- 
ized wave in a plasma has been reported p^ . However, 
angular momentum considerations are not necessary and 
don't seem to be helpful for the derivation of the IFE. 



VI. CONCLUSION 

The ongoing research efforts to find faster and more 
precise methods to store information in increasingly 
miniaturized magnetic particles has lead to several suc- 
cessful collaborations between scientists specialized in 
different domains. This is especially true for modern 
experimental investigations, e.g., on ultrafast magneti- 
zation process in nanostructures, where usually several 
different techniques are involved. The recently demon- 
strated manipulation of spins by means of a circularly po- 
larized laser beam is a further advance in this field, which 



4 



requires an expansion of the knowledge about magnetiza- 
tion processes driven by electromagnetic waves. Similar 
to the positive synergy effects that have been obtained 
in experiments, the theory can also benefit significantly 
by combining the theory of magnetism with theoretical 
plasma physics, in which the interaction of electrons with 
radiation is of central importance. 

Compared to the relatively well known phenomenolog- 
ical treatment of the IFE, the microscopic derivation of 
the IFE presented in this paper gives a clearer insight 
into the processes leading to the magnetization of the 
sample. A further important aspect is that eq. 1)24(1 is 
free of material constants and thus allows for quantitave 
predictions of the IFE. In particular, a strong depen- 
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